Keywords

1 Introduction

Ultra-precise inter-satellite ranging as performed for more than 15 years by the GRACE mission has been established as the state-of-the-art technique to globally observe mass variations in the system Earth from space (Tapley et al. 2019). Continued meanwhile by its Follow-On mission (GRACE-FO, Flechtner et al. 2013), a growing number of institutions is processing the GRACE/GRACE-FO Level-1B instrument data to derive mass variations on a monthly basis (Level-2 products). Although each new release of monthly gravity fields represents a significant improvement with respect to earlier releases, the solutions of different institutions usually differ considerably in terms of noise (Jean et al. 2018) and sometimes also in terms of signal (Meyer et al. 2015). In the frame of the European Gravity Service for Improved Emergency Management (EGSIEM) initiative, a prototype of a scientific combination service has been set up to demonstrate that improved solutions may be derived by combining individual solutions which are based on different approaches but agreed-upon processing standards (Jäggi et al. 2019). The Combination Service of Time-variable Gravity Fields (COST-G) continues the activities of the scientific combination prototype service of the EGSIEM initiative to realize a long-awaited standardization of gravity-derived mass transport products under the umbrella of the International Association of Geodesy (IAG). Established at the 2019 General Assembly of the International Union of Geodesy and Geophysics (IUGG) as a new Product Center of IAG’s International Gravity Field Service (IGFS) for time-variable gravity fields, COST-G will operationally provide consolidated monthly global gravity models with improved quality, robustness, and reliability both in terms of spherical harmonic (SH) coefficients and thereof derived grids of surface mass changes by combining solutions or normal equations (NEQs) from COST-G analysis centers (ACs). The COST-G ACs are adopting different analysis methods but apply agreed-upon consistent processing standardsFootnote 1 to deliver time-variable gravity field models, e.g. from GRACE or GRACE-FO low-low satellite-to-satellite tracking (ll-SST) or from non-dedicated data such as GPS high-low satellite-to-satellite tracking (hl-SST) or Satellite Laser Ranging (SLR). In addition COST-G makes use of existing and publicly available solutions or NEQs of Partner Analysis Centers (PCs), who are directly linked with the GRACE and GRACE-Follow On project. PCs are producing quality controlled products following their own project requirements which may not necessarily be in compliance with the COST-G standards.

The article is structured as follows: Sect. 2 provides an overview of the COST-G organization, Sect. 3 describes the COST-G workflow by discussing the first release of COST-G GRACE Level-2 products and Sect. 4 concludes the article with a summary and future perspectives.

2 COST-G Organization

COST-G is organized in close analogy to other IAG services (Drewes et al. 2016). It is steered by a Directing Board, which sets the objectives and the scientific and operational goals. The COST-G objectives to derive time-variable gravity fields with improved quality, robustness, and reliability are accomplished by the following components.

2.1 Central Bureau (CB)

The Central Bureau is responsible for all operational activities of the Service. The Central Bureau coordinates COST-G activities, facilitates communications and maintains documentations. The CB is currently located at the Astronomical Institute of the University of Bern (AIUB).

2.2 Analysis Centers (ACs)

The COST-G Analysis Centers produce time-variable gravity field solutions according to the specifications defined by the COST-G Processing Standards defined by the COST-G Directing Board. For the analysis of GRACE/GRACE-FO data the current ACs are, in alphabetical order:

  • Astronomical Institute, University of Bern (AIUB)

  • Centre National d’Etudes Spatiales (CNES/GRGS)

  • German Research Centre for Geosciences (GFZ)

  • Institute of Geodesy, Graz University of Technology (IfG/ITSG)

The list of ACs may differ for the processing of non-dedicated satellite data to derive alternative monthly solutions of the Earth’s time-variable gravity field, e.g. from Swarm hl-SST data in the frame of an ESA initiative (Teixeira da Encarnação et al. 2019).

The ACs send their solutions to the Analysis Center Coordinator for combination together with a summary describing their processing strategy. Depending on the availability of new or improved AC contributions, new combined solutions are released. GRACE contributions need to cover at least the time period 2003 to mid 2016 to be included in the combination. Shorter periods may be considered for testing. Corresponding rules for GRACE-FO, hl-SST, SLR contributions will be defined when a significantly large number of ACs is available and has consolidated their processing strategies.

2.3 Partner Analysis Centers (PCs)

COST-G will in addition make use of existing and publicly available solutions or NEQs of other processing centers, denoted as Partner Analysis Centers (PC). Currently, these are primarily centers that are part of the GRACE and GRACE-FO project, e.g. the Center for Space Research (CSR) of the University of Texas at Austin or NASA’s Jet Propulsion Laboratory (JPL), who are producing quality controlled products following their own project requirements. COST-G retains the right to exclude solutions if they either deviate from the COST-G Processing Standards or do not comply with decisions of the COST-G Directing Board.

2.4 Analysis Center Coordinator (ACC)

The Analysis Center Coordinator, currently located at AIUB, first compares the individual gravity fields for quality control. From the comparison of the solutions passing the quality control the ACC defines empirical weights for the individual contributions. Eventually, he combines the accepted solutions to generate a combined gravity field using the underlying normal equations of the individual ACs (Meyer et al. 2019). If normal equations are not available, combinations may also be performed on the solution level (Jean et al. 2018). The resulting COST-G solutions (Level-2 products) are published at the International Center for Global Earth Models (ICGEM, Ince et al. 2019).

2.5 Validation Center (VC)

Validation of the COST-G products happens at the Validation Center, which is currently operated by the COST-G ACs at CNES and GFZ. This involves the evaluation of the noise of the solutions over dedicated areas of low variability, as well as the evaluation of the solution quality through comparison with external data from altimetry. Eventually the COST-G products are also validated through fits of tracking data of Low Earth Orbiters.

2.6 Product Evaluation Group (PEG)

External expert users not associated with COST-G ACs are forming the COST-G Product Evaluation Group. They assess the COST-G products for studying mass variations related to terrestrial water storage over non-glaciated regions, bottom pressure variations in the oceans and ice mass changes in Antarctica and in Greenland.

2.7 Level-3 Product Center (L3C)

Various corrections and reductions have to be applied to the Level-2 products, resulting in post-processed SH coefficients (Level-2B products), before user-friendly Level-3 products are generated by the Level-3 Product Center, currently located at GFZ. The COST-G Level-3 products are visualized and described at GFZ’s Gravity Information Service (GravIS, http://gravis.gfz-potsdam.de) and will be made available by GFZ’s Information System and Data Center (ISDC, https://isdc.gfz-potsdam.de). Additionally, the products can also be visualized at the COST-G Plotter (http://cost-g.org/).

3 COST-G Workflow

The workflow of COST-G is illustrated by Fig. 1. It consists of harmonization and quality control of the individual input solution of the different ACs and PCs, and the combination and validation of the resulting COST-G solution. All aspects of this process are explained and illustrated in more detail in the subsequent subsections by discussing the results of the first release of combined GRACE Level-2 products (GRAC_COSTG_BF01_0100, subsequently addressed as COST-G RL01).

Fig. 1
figure 1

COST-G workflow

3.1 Harmonization of Input Solutions

The individual input solutions of the different ACs may differ in the underlying constants (Earth’s gravity constant, Earth’s equatorial radius), tide system, and mean pole convention (see Eq. 21 in Wahr et al. (2015)). Individual solutions may thus need to be re-scaled and individual coefficients (C20, C21, S21) may need to be further corrected to be compliant with the tide system and mean pole convention used by COST-G. Additional background models, such as the atmosphere and ocean dealiasing product (AOD), may also differ between groups and are consolidated before the combination to ensure a consistent signal definition.

3.2 Quality Control of Input Solutions

For quality control of the COST-G products, the signal and noise content of the harmonized AC contributions are first compared. The signal content is analyzed by computing the amplitude of seasonal variations of equivalent water height (EWH) for a large number of river basins and by computing ice mass trends in Greenland and Antarctica. An assessment of the noise is performed by analyzing anomalies, which represent the monthly variability after subtraction of a deterministic model of secular and seasonal variations. The assessment of the seasonal variations of the input solutions for COST-G RL01 is discussed below, whereas the assessment of ice mass trends will be separately discussed in Sect. 3.4.1. The noise assessment is not reproduced here as it largely corresponds to earlier results presented in Jean et al. (2018) and Meyer et al. (2019).

Figure 2 (top) shows the amplitudes, expressed in mean equivalent water height (MEWH), of the estimated seasonal variations for the 500 largest river basins as derived from all AC and PC solutions used for the COST-G RL01. The underlying river basin masks were taken from the Total Runoff Integrating Pathways (TRIP4Footnote 2) model. The corresponding formal errors of the estimation are shown in Fig. 2 (bottom). The analysis shows that no systematic signal attenuation may be observed for seasonal signals in any of the contributing gravity field time-series. The cluster of outlying larger formal errors visible in Fig. 2 (bottom) is related to regions with small seasonal signals and large non-linear trends as occurring, e.g. for regions with accelerated ice mass loss.

Fig. 2
figure 2

Amplitude of seasonal variations (top) and formal errors of amplitudes (bottom) in MEWH (m) for major river basins using a 400 km Gauss filter

3.3 Combination of Input Solutions

Generally, the planned strategy for future COST-G releases is to provide a combination based on the underlying NEQs of the individual ACs according to the methodology presented in Meyer et al. (2019). As NEQs were not available from all centers, the combination was performed on the solution level for the COST-G RL01 as a weighted combination of the SH coefficients using variance component estimation (VCE) (Jean et al. 2018). The underlying AC und PC solutions are, AIUB RL02 (Meyer et al. 2016), GRGS RL04, GFZ RL06 (Dahle et al. 2019), ITSG-Grace2018 (Kvas et al. 2019), and CSR RL06 (Save et al. 2018).

Apart from AIUB RL02, which is still based on the RL02 of the GRACE Level-1B data, all input solutions are based on the most recent RL03 of the GRACE Level-1B data (GRACE Level 1B JPL Release 3.0 2018). For the GRGS solution the underlying NEQs were inverted by the ACC to obtain a solution without regularization.

Figure 3 shows the relative weights assigned to the individual input solutions as determined by VCE. The weights can be interpreted as quality indicators of the solutions as they are inversely proportional to the noise levels of the individual contributions. The highest weights are usually assigned to the ITSG-Grace2018 solutions due to their very low noise level. This is largely related to the sophisticated empirical noise modeling of the ITSG-Grace2018 solution (Kvas et al. 2019) and in accordance with analyses based on earlier ITSG releases (Jean et al. 2018). Lowest weights are usually assigned to the AIUB RL02 solutions due to the use of (meanwhile) outdated Level-1B data and, in particular, since active thermal control of the GRACE satellites was switched off in April 2011 (Tapley et al. 2015), which would have required adaptions in the accelerometer parametrization as shown in Klinger and Mayer-Gürr (2016).

Fig. 3
figure 3

Weights assigned to input solutions

Figure 4 shows median degree amplitudes of anomalies for the combined solution as well as for the input solutions that contributed to the combination for the entire GRACE time period with respect to a linear and seasonal model without applying any filtering. The analysis reveals that in the spectral domain the main gain of the combination is in the range of degrees 15–45. When truncating all gravity fields at order 29 to exclude the effect of the noisy higher-order SH coefficients, which are usually attenuated in applications by post-processing filters, e.g. Kusche (2007), the gain of the combination may even be seen up to about degree 65. The lower noise of the combined solution may also be confirmed in the spatial domain by analyzing the RMS of EWH anomalies over the oceans (not shown).

Fig. 4
figure 4

Median degree amplitudes of residual SH coefficients. Solid: full spectrum, dashed: limited to order m = 0…29

3.4 Validation of Combined Solution

Internal and external validation is performed to ensure the quality of the COST-G solutions by identifying outlying or systematically deviating input solutions.

3.4.1 Internal Validation

Ice mass loss in the polar regions is of enormous societal relevance (Shepherd et al. 2012). Evaluating GRACE mass change time series for the ice sheets in Antarctica and Greenland, as derived from the individual input solutions and the combined COST-G solution, is thus an essential task of the COST-G product evaluation group to detect potential inconsistencies between the input solutions.

3.4.1.1 Greenland Ice Sheet

GRACE-derived mass change time series are derived from all input solutions that contributed to the combination for the entire GRACE time period and from the combined COST-G solution for all drainage basins of the Greenland Ice Sheet (GISFootnote 3) and for the entire Greenland Ice Sheet by adopting the sensitivity kernel approach (Groh et al. 2019 and references in there).

The estimated linear trends of ice mass change for all drainage basins and for the entire Greenland ice sheet agree very well for the different solutions (not shown). Most notably the COST-G time series is characterized by a very favorable noise behavior. Figure 5 shows the noise levels of the mass change time series available over the entire time span in terms of the scaled standard deviation of the derived noise time series for all Greenland Ice Sheet drainage basins and the entire ice sheet following the method of Groh et al. (2019). For the majority of the basins the COST-G time series shows indeed the lowest noise of all solutions used for the combination.

Fig. 5
figure 5

Noise level of the mass change time series for all GIS basins and the entire GIS (basin no. 9)

3.4.1.2 Antarctic Ice Sheet

A similar analysis is performed for the drainage basins of the Antarctic Ice Sheet (AIS, see footnote 3), selected aggregations, and for the entire Antarctic ice sheet. Whereas most of the results of this analysis confirm the level of agreement found for the Greenland Ice Sheet, the linear trends resulting for different gravity field solutions for some drainage basins deserve special attention. Figure 6 compiles for the solutions available over the entire time span the estimated linear trends and error bars in ice mass change for selected drainage basins, the Antarctic Peninsula, East Antarctica, West Antarctica, and for the entire ice sheet. The displayed error bars account for the formal errors of the estimated trends, as well as for leakage errors and errors in the applied model reductions.

Fig. 6
figure 6

Linear trend and error bars in ice mass for drainage basins of the AIS (6, 13, 17, 22), selected aggregations (29: Antarctic Peninsula, 30: East Antarctica, 31: West Antarctica) and the entire ice sheet (32)

Figure 6 shows that the trends resulting from all solutions agree fairly well for West Antarctica. It also reveals, however, that discrepancies are occurring for East Antarctica. Whereas the trends resulting from the CSR RL06 and the ITSG-Grace2018 solutions are in almost perfect agreement, the GFZ RL06 solution suggests a different trend for this entire region. However, all trend estimates agree within their error estimates, since these are clearly dominated by the error of the reduced glacial isostatic adjustment model. As a consequence of the weighting scheme used for the combined solution, the trend of the COST-G solution is in-between the two concurring trends. Due to the larger weights assigned to the CSR / ITSG solution (cf. Fig. 3), it is closer to these solutions. Figure 6 also shows that the different trends over East Antarctica influence the trends resulting for the entire ice sheet. Trend analyses based on shorter time periods common to all solutions show again different trends for the AIUB RL02 and GRGS RL04 solutions (not shown).

3.4.2 External Validation

Currently two methods are realized within COST-G to assess the quality of the solutions by independent data sets. More external validations may follow in the future.

3.4.2.1 Comparison to Altimetry

Currently two test areas (Caspian Sea, Black Sea) are selected within COST-G for an independent signal assessment. The time series of the time-variable gravity field solutions are compared with the time series of altimetric heights, derived from Hydroweb for the Caspian Sea and AVISO+ for the Black Sea. One bias and one scale factor are adjusted to perform the comparison.

Table 1 exemplarily shows the correlation coefficients over the Caspian Sea when filtering the gravity field solutions with different DDK filters (Kusche 2007). It can be seen that the COST-G solution presents in this metric a slight improvement with respect to the input solutions.

Table 1 Correlation with altimetry over the Caspian Sea
3.4.2.2 Orbital Fits

The long to medium wavelength accuracy of gravity field models can be evaluated through dynamic orbit computations as commonly done for the evaluation of static gravity fields, e.g. Gruber et al. (2011). In the frame of COST-G the same concept is adopted to the time-variable gravity field solutions from the individual ACs and thereof derived combined solution. In order to not suffer from large omission errors of the time-variable GRACE solutions, which are generally only available up to degree and order 90, all solutions are filled up to degree and order 240 with the SH coefficients of the GOCE static model DIR-6 (Förste et al. 2019). For the COST-G RL01 dynamic GOCE orbits with an arc length of 1.25 days were fitted to the GOCE kinematic precise science orbits (Bock et al. 2014) serving as pseudo-observations. The gravitational forces were modeled according to the different gravity field models under consideration, whereas non-gravitational accelerations were described by the high-quality GOCE accelerometer data. For each arc three common mode acceleration biases are estimated in addition to the initial state vector. The scale factors of the common mode accelerations were fixed to one (Gruber et al. 2011).

Table 2 shows the RMS of orbit fits for the different test cases, derived as mean values from the 3D residuals from the 30/31 individual arcs under consideration for each month. It can be seen that the COST-G RL01 solution is also assessed in this metric of very good quality, but there is potential for further improvements as the solution is not yet best for each of the tested months. Note that the orbit test does not primarily validate C20 of the individual solutions. Only marginal differences would result if C20 is replaced in all solutions with one and the same value (not shown).

Table 2 RMS of dynamic GOCE orbit fits (cm) for the different gravity fields

4 Conclusions

COST-G has officially started its operational service at the IUGG’s 2019 General Assembly as a new Product Center of IAG’s International Gravity Field Service (IGFS). It continues the scientific combination prototype service that has been established within the EGSIEM initiative to realize a long-awaited standardization of gravity-derived mass transport products under the umbrella of the IAG. COST-G recognizes and emphasizes the existence and acknowledges the contribution of every individual analysis center and partner analysis center to this community effort. Their participation is a crucial and mandatory prerequisite to the consolidation of monthly global gravity fields within COST-G.

At the occasion of the 2019 IUGG General Assembly COST-G has provided a first release of combined GRACE monthly solutions covering the entire GRACE time period between April 2002 and June 2017 by combining the solutions of five contributing centers. In addition COST-G also provides combined Swarm monthly solutions as an operational product in the frame of an ESA initiative (Teixeira da Encarnação et al. 2019). Depending on the interest of the scientific community, further products may be established in the future, e.g. combined monthly solutions derived from SLR satellites. For GRACE and GRACE-FO it is planned to base future COST-G releases whenever possible on the combination of the normal equations of the underlying input solutions.

COST-G has set up a workflow that allows for a rigorous evaluation of the products by both internal and external means. The experience from the first release of combined monthly GRACE solutions has underlined that such a rigorous evaluation is of key importance to ensure the quality of the COST-G products and will therefore be further extended in the future. Although the COST-G RL01 solutions show an excellent behavior in terms of noise, as it is demonstrated by the various internal and external quality metrics, a discrepancy in the signal has been revealed by the COST-G product evaluation group for the different input solutions when regarding ice mass trends over East Antarctica. As this also affects the trends of the COST-G RL01 for that region, the root cause is currently being further investigated. The COST-G workflow also includes the generation of Level-2B products to derive thereof user-friendly Level-3 products that will enable hydrologists, glaciologists, oceanographers, geodesists and geophysicists to fully profit from one combined, and consolidated monthly GRACE gravity product. Updates will be announced at http://cost-g.org.